The magnetic field dependence of the electronic specific heat of Y o s Cao.2Ba2CusOe +x 
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We have measured the electronic specific heat of Yo.gCao.2Ba2Cu307-s using a high-resolution 
differential technique from liquid helium temperature to room temperature in an applied magnetic 
field up to 13T. The field dependence of the electronic specific heat at low temperatures in the 
superconducting state behaves differently in the overdoped and underdoped regimes, varing as \fK 
in the overdoped regime but as H in the underdoped regime. An entropy loss is observed in the 
normal state in optimal and underdoped samples, which can be attributed to the opening of a normal 
state psuedogap. From the temperature and field dependences of the free energy, the temperature 
dependence of the upper critical field H C 2 is determined. For the overdoped sample (x = 0.79), we 
find H C 2 to have a negative curvature. 



I. INTRODUCTION 

Measurements of electronic specific heat have played 
an important role in clarifying the nature of supercon- 
ducting pairing symmetry , normal state pseudogap 
, quantum criticality |)]| , and many other fundamen- 
tal properties of high-T c superconductors Jh|[ii| ■ While 
the mechanism of high-temperature superconductivity is 
not understood, there is an emerging consensus that the 
superconducting state itself is not particularly exotic, in 
the sense that it consists of a BCS-like pair state with 
well-defined quasiparticle excitations. Specific heat ex- 
periments are insensitive to the phase of the order param- 
eter, but can provide bulk information on the behavior of 
low energy density of states. It has been confirmed that 
at very low temperatures below T c the zero field specific 
heat varies as T 2 and the field dependent specific heat 
scales as y/~H, fl],|],[l consistent with the d-wave pair- 
ing scenario of high-T c superconductivity. Revaz et.al |jj 
also showed that the field dependent term obeys a scal- 
ing function of T / \f~B proposed by Volovik jL2[ . The loss 
of the entropy deduced from the specific heat data well 
above T c in underdoped cuprates has revealed much im- 
portant information on the normal state pseudogap |^| || . 

In this paper we report our experimental data of the 
electronic specific heat of Yo.sCio^-Ba^CuaOg+z mea- 
sured from ~ 7K to room temperatures in an external 
magnetic field up to 13T. We choose Calcium doped 
Y Ba2CuzOQ+ x (YBCO) samples because the substitu- 
tion of Y 3+ by Ca 2+ adds holes to the material. This 
allows physical properties of YBCO-type cuprates to be 
investigated in both over- and under-doped regimes. 

The data reported in this paper were taken on poly- 



crystalline Yq&C 'ao. 2 BCI2CU3O6+X samples of about lg 
in weight each, prepared by conventional solid state reac- 
tion. The oxygen content was controlled by the annealing 
temperature and determined from the weight loss after 
annealing, with an uncertainty 0.01 in x. We measured 
the specific heat using a differential technique in which 
the difference in specific heat between the sample under 
investigation and a related reference sample is measured 
with a resolution of about 1 : 10 4 (i.e. 1% of the elec- 
tronic contribution) [p"3f . The majority of the phonon 
contribution is backed off by the reference sample and a 
simple and reliable correction for residual differences in 
phonon terms allows us to separate accurately the elec- 
tronic specific heat C e i from the total ||. In this paper 
we are mostly concerned with the change of the specific 
heat between a finite field and a zero field determined 
in separate runs. Whilst this quantity is independent of 
the phonon specific heat, the inherent reproducibility of 
the differential technique coupled with the elimination of 
the bulk of the phonon term greatly increases the relia- 
bility with which small field dependent changes can be 
determined. 

The rest of the paper is arranged as follows. In Sec. 
||, our experimental results are presented and discussed. 
Emphasis is given to the scaling behavior of the specific 
heat in a finite field and the temperature dependence of 
the upper critical field. Sec. Ill gives a brief summary. 



II. RESULTS AND DISCUSSIONS 

We have measured seven 20% Ca substituted YBCO 
samples. Zero field measurements on this system have 
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been discussed previously 0. Fig.[IJ shows the critical 
transition temperature T c for all the samples. The op- 
timal doping is located at x = 0.7 and T c = %AK. The 
condensation energy at zero temperature f (0), rescaled 
by T 2 , is also shown in this figure. The value of {7(0), 
as will be discussed later, is determined from the dif- 
ference of the free energy density between the normal 
and superconducting states. C/(0)/T c 2 is nearly doping 
independent when x > 0.79 but drops rapidly with de- 
creasing x when x < 0.79. U(0)/T 2 at the optimal dop- 
ing x = 0.7 (T c = 84K) is about half of the value at 
x = 0.79 (T c = 80K). This rapid decrease in {7(0) is due 
to the opening of the normal state pseudogap at a planar 
hole density p crit ~ 0.19holes/CuO2 The abrupt 

change in the doping dependence of the condensation en- 
ergy around x — 0.79 might be a manifestation of quan- 
tum critical phenomena. Similar changes have also been 
found in the doping dependence of the supcrfluid stiffness 
as well as many other physical quantities [j9 15|. 
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FIG. 1. Superconducting transition temperature T c and 
(7(0)/T c 2 for Y . & Ca . 2 Ba 2 Cu 3 O 6+x , where U(0) is the con- 
densation energy at zero temperature. The optimal doping is 
located at x = 0.7 with T c = 8AK. {7(Q)/T C 2 is nearly dop- 
ing independent when x > 0.79 but decreases with decreasing 
doping when x < 0.79. 



A. Specific heat and the scaling behavior 

The specific heat is an important quantity in character- 
izing low-energy excitations (see for example ref. pfl| , pT[ |). 
It provides a direct measure of thermal excitations (spec- 
tral weight) in the electronic spectrum. In a Fermi liq- 
uid system, the temperature dependence of the electronic 
specific heat coefficient, ^(H,T) = C e i(H,T)/T at low 
temperatures reflects the energy dependence of the den- 
sity of states of low-lying excitations. 

Fig. |^ shows the change of the specific heat coefficient 
Aj(H, T) = j(H, T) - 7 (0, T) in four fields H = 2T, AT, 
9T and 13T for Yo^Cao,2Ba2Cu30ejg. Since no correc- 
tions from the lattice contribution to the raw data of A7 
are involved, there is no ambiguity in this quantity. In 




the normal state (T > T c + 1QK) , 7 is unchanged by 
the applied fields within an experimental uncertainty of 
~ 0.01m J I gat. K 2 . At zero field (see the inset of Fig||, 
7 shows a sharp peak around T c . This peak is gradually 
suppressed by increasing field. 
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FIG. 2. Changes of the electronic specific heat co- 
efficient Ay(H, 0) at four different magnetic fields for 
lo.8Cao.2Ba2C1i3O6.79. H = 2T, AT, 9T, and 13T from bot- 
tom to top. 

In type-II superconductors in an applied magnetic 
field, there are two types of low-lying excitations: vor- 
tex core bound states and extended quasiparticle states 
which evolve smoothly into the bulk zero-field quasipar- 
ticle states at large distances from the vortex |l2] ]. In a 
conventional s-wave superconductor, the low energy den- 
sity of states is dominated by the core bound states, since 
the extended states are fully gapped and the core level 
spacing is about a factor A/sf smaller than the quasi- 
particle excitation gap, where A is the energy gap and ep 
is the Fermi energy. Since the number of vortices scales 
linearly with H , the low energy density of states is also 
proportional to H . This means that at low temperatures 
and when H <C H C 2, the specific heat scales linearly with 
H/H C2 , i.e. C el {H)~H. 

However, in the mixed state of a superconductor with 
line nodes, Volovik |L2| pointed out that the low energy 
density of states is dominated by contributions from ex- 
tended quasiparticle states rather than the core bound 
states. It was shown that the core bound states con- 
tribute less to the entropy than the extended quasiparti- 
cles by a factor that diverges logarithmically with the 
inter-vortex separation. This is due to the effect of 
Doppler shift caused by supercurrents circulating around 
a vortex. Since the inter-vortex distance is proportional 
to i7 -1 / 2 and the number of vortices is proportional to H, 
it follows that the density of states at the Fermi level is 
proportional to y/~H and the specific heat varies as T\f~H 
in the limit T — > 0. Furthermore, from the \f~H. scaling 
behavior of the excitation spectrum, Simon and Lee JlfJ] 
showed that low temperature thermal and transport co- 
efficients of a d-wave superconductor should scale with 
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<JH/T. In particular, they predicted that the specific 
heat has a simple scaling form 



c el = tVhf c (t/Vh 



(1) 



where Fq is a undetermined scaling function. This scal- 
ing prediction holds for clean d-wave superconductors 
and for energy scales small compared to the maximum 
gap scale. In the dirty limit, Kubert and Hirschfeld p"^ ] 
showed that the low-lying excitations are still dominated 
by the extended states, but the specific heat behaves as 
HlnH at low fields. The \[HT term of the specific heat 
was first identified by Moler et al R] for YBCO single 
crystals and later by other groups [fpjj. The \[H JT scal- 
ing behavior has also been reported in the specific heat 
measurements of YBCO 
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FIG. 3. Low temperature behaviours of Ay(H,T) at five 
different fields for Yq,%C ao,2Ba2CuzO§.7$. The solid lines are 
linear fits to the data. 

For Ca doped YBCO, we find that the field dependent 
A7 varies linearly with temperature at low temperatures 
(Fig. ||) and can be fitted by the equation 



A-y(H, T) w A 1 (H 1 0) + a(H)T 



(2) 



within experimental errors. The first term is just the co- 
efficient of the INITIAL linear T term of C e i and is pre- 
dicted to vary as \[H in an ideal d-wave superconductor. 
For overdoped YQ^Cao.2Ba2Cu^,OQ+ x (the upper panel 
of Fig. g), we indeed find that Aj(H,0) ~ y/H, with 
slopes comparable with that found for optimally doped 
single crystal YBCO and the theoretical predication 
(P~2y . However, for underdoped Yo,sCao.2Ba2Cu30e+ x 
(the lower panel of Fig. ||), we find that A^{H, 0) varies 
linearly with H. This is different from the \J~H scaling 
behavior of a clean d-wave superconductor and also from 
the HlnH scaling behavior predicted for a dirty d-wave 
superconductor. We do not know if this violation of the 
\f~H~ scaling law of C e i in these underdoped samples is 
due to the effect of normal state pseudogap or just be- 
cause the temperatures we have measured are not low 



enough to probe the characteristic low-energy feature of 
a d-wave superconductor. Further investigation on the 
field dependence of the specific heat of underdoped high- 
T c materials at low temperatures is desired. 
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FIG. 4. The specific heat coefficient A*y(H, 0) as a 
function of \[H for the overdoped Yo.sCao.2Ba2Cu30e+ x 
(upper panel) and a function of H for the underdoped 
YQ.sCao.2Ba2Cu2,0§+ x (lower panel). 

The field dependence of ct(H) is more complicated and 
less accurately determined than Aj(H). However, for the 
x = 0.79 sample, we find that a(H) varies almost linearly 
with H (Fig. |^). This linear H dependence of ct(H) is 
not what one may expect from the scaling relation (|lj). 
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FIG. 5. a{H) versus H for Y , g Cao.2Ba2Cu 3 Ofj.7 
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B. Entropy 



The electronic entropy S is one of the parameters sen- 
sitive to both superconducting condensation and normal 
states pseudo-gap and can be determined from the tem- 
perature integration of thespecific heat coefficient 



S 



l{T)dT 



(3) 



Our data of S for 20% Ca doped YBCO at zero field 
are shown in Fig. ^[ There is an abrupt change in the 
slope of S at T c due to the specific heat jump. For the 
overdoped samples (x=0.93, 0.87 and 0.79), the normal 
state entropy extrapolates to zero at zero temperature. 
However, for the optimal and underdoped samples (x= 
0.7, 0.64, 0.56 and 0.48) the normal state entropy ex- 
trapolates to a negative value at zero temperature. This 
means that some of the entropy is lost in the normal 
state in these materials. This loss of entropy has been 
attributed to the opening of the normal state pseudo- 
gap. Similar results have been observed in other cuprate 
systems. |6|| 
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FIG. 6. Electronic entropy versus temperature for 
Yo,gCao,2Ba2Cu^06+ x at zero field, x for the curves from 
top to bottom are 0.93, 0.87, 0.79, 0.70, 0.64, 0.56, and 
0.48, respectively. The dashed line is for the optimal dop- 
ing (x=0.70). The entropy loss is observed in normal state in 
optimal as well as underdoped samples. 

Fig. shows the enhancement of the entropy at a 
13T magnetic field for different oxygen contents. The su- 
perconducting condensation is suppressed by an external 
field below T c . This causes the increase of the entropy. 
In the normal state, the entropy is almost unchanged by 
the external field. This suggests that the normal state 
pseudogap is unaffected by an external field, consistent 
with the NMR data for underdoped YBCO El. 
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FIG. 7. Difference of the entropy between zero field and 
H — 13T for Yo.gCao.2Ba,2Cu306+x. The value of x is given 
besides each curve. The external field affects strongly the 
condensation energy but has almost no effect on the normal 
state entropy. 
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FIG. 8. (a) AF N . S (T) = F N (0,T) - F S (Q,T) versus 
temperature for Yo.sCao.2Ba2Cu30e+ x . (b) AFs(H,T) = 
F S (H, T) - F s (0, T) at six different fields for x = 0.79. 



C. Free energy 

From the entropy, one can further find out the free 
energy density. In particular, by integrating out the en- 
tropy difference between the normal and superconducting 
states at zero field, the change of the free energy 
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AFjv, s (T) = F N (0,T) - F S (P,T) 



8tt 



(4) 



can be obtained, where H C (T) is the thermodynamic 
field. At zero temperature AFn,s(0) is the condensa- 
tion energy U(0). Our results of AFn.s{T) are shown in 
Fig. ||(a) for different oxygen contents. The behavior of 
this quantity is different for overdoped and underdoped 
samples. With decreasing temperatures, AF N ^ S (T) in- 
creases more rapidly in the overdoped samples than in 
the underdoped ones. 

The change of the free energy density in a field in the 
superconducting state AF S {H,T) = F s (H,T)-F s (0,T) 
can be similarly obtained. Fig. ||(b) shows AFs(H,T) 
for the x — 0.79 sample. The largest AF$ at 13T and 
OK is about 220mJ/gat, much smaller than £7(0) ~ 
1550m J/g. at. This indicates that the upper critical field 
H C 2 at zero temperature is much larger than our maxi- 
mum field of 13T. 



the free energy AF$ in the superconducting state can be 
deduced from the reversible magnetization M: 



AF S (H,T) 



H 



MdH. 



(5) 



The London formula of the magnetization in the mixed 
states is approximately given by [E7| : 



M = - 



87rm/ 



In 



bH, 



<::2 



H 



(6) 



where k = A/£ is the Ginsberg-Laudau parameter, A is 
the penetration depth, and & is a constant. The lower 
and upper critical fields H c i and 7J c2 are related to the 
thermodynamic field by the formula H c iH C 2 = H^Iuk. 
Substituting (^J) into (j^) and carrying out the integration, 
we find that 



AF S (H,T) 



HAF N , S (T) 
H c2 (T) 



, bH ^ i 1 



(7) 



D. The upper critical field 

In recent years, there has been intensive investigation 
on the temperature dependence of the upper critical field 
H C 2 in high-T c superconductors. MacKenzie et al. JToj ] 
measured the field dependence of the in-plane resistiv- 
ity of overdoped T^i^CuOg and found that H C 2(T) 
shows a positive curvature at low temperatures. Simi- 
lar temperature dependence of H C 2 has also been found 
in the transport measurements of other overdoped ma- 
terials p0|| . To understand this puzzling upward curva- 
ture of H c2 (T) , a number of theories have been proposed 
pl| invoking thermal and quantal fluctuations. Re- 
cently, Wen et al measured temperature dependence of 
the magnetic moment of overdoped Bi2Sr2- x La x CuOi 
and La,2-xSr x Cu04 L and found that there are two tran- 
sitions in the M(T) curve in an applied field J2^,||F|]. The 
field dependence of the lower transition (or crossover) 
temperature indeed exhibits an upturn at low tempera- 
tures, similar to the upper critical field found in the trans- 
port measurements. However, the field dependence of the 
upper transition temperature behaves very differently. It 
shows a negative curvature, same as for the upper critical 
field in conventional superconductors. They argued that 
the transition at higher temperature corresponds to the 
true upper critical field H C 2 and the upward curvature 
of the critical field at lower temperature is in fact due to 
the formation of a bulk superconducting phase coherence 
through Josephson coupling between some superconduct- 
ing clusters with a transition temperature higher than the 
bulk transition temperature 21 1. 

We can estimate the upper critical field H c2 from the 
free energy data obtained above if we assume that the 
magnetisation obeys the London expression (see for ex- 
ample Ref pi). In the mixed state, the difference of 
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FIG. 9. AFs(H,T) versus H at different temperatures for 
Yo.sCao.2Ba2Cu3Oe.7g- The dot and dash lines are the fits 
using Eq.(7) 

Fig. H shows AF S ) as a function of H at different tem- 
peratures for the x = 0.79 sample. The the dash and 
dot lines are the fits of AF$ using Eq.(^) with the value 
of Fn.s(T) shown in fig.||. The values of H C 2 for this 
material can be obtained from the fits and the result is 
shown in Fig. 10 for the x = 0.79 sample (similar re- 
sults have been obtained for other samples) . Same as for 
conventional s-wave superconductors, we find that H C 2 
has a negative curvature. Close to T c , H C 2 varies lin- 
early with temperature and —dH C 2/dT ~ bT/K. This 
agrees with the electronic specific heat measurement for 
overdoped Tl2Ba 2 C'uOe |2£|, where H c2 is determined 
by the changes of T c in magnetic fields. The magni- 
tude of H C 2 is much larger than that determined from 
the transport measurement, in agreement with the mag- 
netization measurement data of Wen et al |25U26f| for 
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overdoped BiiSr-2.- x La x CuO '4 and La.2- x Sr x CuO±, but 
in disagreement with the transport measurement data 
for overdopedTZ 2 Ba 2 Cu0 6 |l9| and Bi 2 Sr 2 CuO y |o). 
The difference between the thermodynamic and trans- 
port measurements has been attributed by Wen et al to 
a macroscopic phase separation of electrons in overdoped 
materials. Another alternative is that transport mea- 
surements are sensitive to flux-flow resistance at tem- 
peratures much lower than the thermodynamic critical 
temperature. 
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FIG. 10. The upper critical field H c2 for 

lo.8Cao.2-Ba2C1t3O6.79. 



III. SUMMARY 

In summary, we have measured the electronic specific 
heat of Yo,sCao.2Ba2Cu30e+ x using a differential tech- 
nique from IK to room temperature in a magnetic field 
up to 13T. We have discussed the physical properties 
of the specific heat, the entropy, the free energy density 
and the upper critical field. For optimal and overdoped 
samples, we found that the field dependent specific heat 
coefficient scales with VH at low temperatures, in agree- 
ment with other experimental measurements. However, 
for underdoped samples, we found that the low tempera- 
ture specific heat coefficient scales linearly with H, incon- 
sistent with the \fli scaling law as expected for clean d- 
wave superconductors. Further investigation on the field 
dependence of low temperature specific heat in clean un- 
derdoped high-T c superconductors should be carried out 
to resolve this discrepancy. The upper critical field for 
the overdoped sample (x — 0.79) shows a downward cur- 
vature. This agrees with the results of the magnetiza- 
tion measurements but differs from those of the transport 
measurements. It suggests that the overdoped materials 
might be very inhomogeneous in an external field and 
the upward curvature of H C 2 observed in the resistivity 
measurement is probably due to the Josephson coupling 
between the superconducting clusters produced by elec- 
tronic phase separation. 
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